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Abstract 

Microgauss magnetic fields are observed in all galaxies at low and high redshifts. The origin 
of these intense magnetic fields is a challenging question in astrophysics. We show here that the 
natural plasma fluctuations in the primordial universe (assumed to be random), predicted by the 
Fluctuation-Dissipation-Theorem, predicts ~ 0.034 fj,G fields over ~ 0.3 kpc regions in galaxies. 
If the dipole magnetic fields predicted by the Fluctuation-Dissipation-Theorem are not completely 
random, microgauss fields over regions > 0.34 kpc are easily obtained. The model is thus a strong 
candidate for resolving the problem of the origin of magnetic fields in < 10 9 years in high redshift 
galaxies. 
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I. INTRODUCTION 



The origin of large-scale cosmic magnetic fields in galaxies and protogalaxies remains a 
challenging problem in astrophysics jl-4]. There have been many attempts to explain the 
origin of cosmic magnetic fields. One of the first popular astrophysical theories to create 
seed fields was the Biermann mechanism j^J. It has been suggested that this mechanism acts 
in diverse astrophysical systems, such as large scale structure formation js-8], cosmological 
ionizing fronts 9j, star formation and supernova explosions 10|, Ryu et al. 12] made 
simulations showing that cosmological shocks can create average magnetic fields of a few 
fiG inside cluster/groups, ~ 0.1 [iG around clusters/groups, and ~ 10 nG in filaments. 



Medvedev et al. 



13j | showed that magnetic fields can be produced by collisionless shocks in 



galaxy clusters and in the intercluster medium (ICM) during large scale structure formation. 
Arshakian et al. [l^| studied the evolution of mag netic fields in galaxies coupled with 



15] investigated second-order couplings 



hierarchical structure formation. Ichiki et al. 
between photons and electrons as a possible origin of magnetic fields on cosmological scales 
before the epoch of recombination. The creation of early magnetic fields generated by 



cosmological perturbations have also been investigated 16l-ll9|. 



In our galaxy, the magnetic field is coherent over kpc scales with alternating directions 



in the arm and inter-arm regions (e.g.,Kronberg 



expected for magnetic fields of primordial origin 22] . 



20], Han 21]). Such alternations are 



Various observations put upper limits on the intensity of a homogeneous primordial mag- 
netic field. Observations of the small-scale cosmic microwave background (CMBj anisotropy 



23|. Reion- 



yield an upper comoving limit of 4.7 nG for a homogeneous primordial field 
ization of the Universe puts upper limits of 0.7 — 3 nG for a homogeneous primordial field, 
depending on the assumptions of the stellar population that is responsible for reionizing 



the Universe 



24] . Another upper limit for a homogenous primordial magnetic field is the 



26] . Thus, if we are investigating the collapse 



magnetic Jeans mass ~ 10 10 M Q (B/3nG) 3 (25|, 
of a ~ 1O 7 M protogalaxy, the homogeneous primordial magnetic field must be < 0.3 nG 
in order for collapse to occur. 

Galactic magnetic fields have been suggested to have evolved in three main stages. In the 
first stage, seed fields were embedded in the protogalaxy. They may have had a primordial 
origin, as suggested in this paper. Another possibility is that the seed fields could have been 
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injected into the protogalaxies by AGN jets, radio lobes, supernovas, or a combination of 
the above. Still another possibility is that the seed fields may have been created by the 
Biermann battery during the formation of the protogalaxy. In the second stage, the seed 
fields were amplified by compression, shearing flows, turbulent flows, magneto-rotational 
instabilities, dynamos or by a combination of the above. In the last stage magnetic fields 



were ordered by a large scale dynamo [27j. 

Ryu et al. 12| investigated the amplification of magnetic fields due to turbulent vorticity 
created at cosmological shocks during the formation of large scale structures. A given 
vor ticity u can be characterized by a characteristic velocity V c over a characteristic distance 



Ryu et al. 



found that u typically is 

U-1-3X10-V 1 , (1) 



which corresponds to 10-30 turnovers in the age of the universe. They investigated L c > 1 
Mpc hr 1 . We investigate L c ~ 200 kpc h^ 1 in protogalaxies for a similar vorticity. 

We show that a seed field 0.003 nG over a comoving 2 kpc region at z ~ 10, predicted 
by the Fluctuation-Dissipation Theorem jsj, amplified by the small scale dynamo is a good 
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Subramanian 



candidate for the origin of magnetic fields in galaxies. K. Subramanian 
29 1 and Brandenburg & Subramanian 3o| derived the non-linear evolution equations for 
the magnetic correlations. We use their formulation for the small scale dynamo and solve 
the nonlinear equations numerically. In §11, we review the creation of magnetic fields due 
to electromagnetic fluctuations in hot dense equilibrium primordial plasmas, as described in 
our previous work j^J. In §111, we discuss the small scale dynamo and in §IV, the important 
parameters of the plasma to be used in the calculations. In §V, we present our results and 
in §VI our conclusions. 

II. CREATION OF MAGNETIC FIELDS DUE TO ELECTROMAGNETIC FLUC- 
TUATIONS IN HOT DENSE PRIMORDIAL PLASMAS IN EQUILIBRIUM 

Thermal electromagnetic fluctuations are present in all plasmas, including those in ther- 
mal equilibrium. The level of the fluctuations is related to the dissipative characteristics 



of the plasma, as described by the Fluctuation-Dissipation Theorem (FDT) 3l| [see also 
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Akhiezer et al. [32|, Dawson 33[, Rostoker et al. 34|, Sitenko 35]]. 

de Souza & Opher ^ studied the evolution of these bubbles as the Universe expanded 
and found that the magnetic fields in the bubbles, created originally at the quark-hadron 
phase transition (QHPT), had a value ~ 9 /iG and a size 0.1 pc at the redshift z ~ 10 (see 
Table 1 of {^j]). Assuming that the fields are randomly oriented, the average magnetic field 
over a region D is B = 9/iG (O.lpc/D) 3 / 2 . The theory thus predicts an average magnetic 
field 0.003 nG over a 2 kpc region at z ~ 10. We assume this seed field and examine its 
amplification in a protogalaxy by the small scale dynamo, discussed in the next section. 



III. SMALL SCALE DYNAMO 

In a partially ionized medium, the magnetic field evolution is governed by the induction 
equation 

(dB/dt) = V x (v. t x B — 77V x B), (2) 

where B is the magnetic field, Vj the velocity of the ionic component of the fluid and 77 is 
the ohmic resistivity. 

Let L c be the coherence scale of the turbulence. Consider a system whose size is > L c 
where the mean field, averaged over any scale, is negligible. We take B to be a homogeneous, 
isotropic, Gaussian random field with a negligible mean average value. For equal time, the 
two point correlation of the magnetic field is 



[B i (x,t)Bl(y,t)) = M i *(r,t) 



(3) 



where 



M 



M, 



N 



28 



S *3 



+ M L 



1 qr 3 



He ijk r 



(4) 



-|30|. Mi(r,t) and (r,t) are the longitudinal and transverse correlation functions, 
respectively, of the magnetic field and H (r, t) is the helical term of the correlations. Since 
V • B = 0, we have = (l/2r) d (r 2 Mi) / (dr) j^ . The induction equation can be 
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converted into evolution equations for Ml and H : 



at r 4 or \ or 

+ G(r)M L (r,t) + 4a N H (r,t) , 



(5) 



and 



dH l_d_ T a d 

a N (r,t)M L (r,t)]] 



r'—[2K N (r,t)H(r,t) 

ar 



(6) 



where 



and 



28 



k n (r, t) = 7] + T LL (0) - T LL (r) + 2 a M L (0, t) 
a N {r,t) = 2 C (0) - 2 C (r) - 4 a if (0, t) , 



CM 



(- 






r 



Id 



(7) 
(8) 

(9) 



29| . Tii(r) and TjvAr(r) are the longitudinal and transverse correlation functions for the 



velocity field. The functions Tjvtv and Tll are then related in the way described by Subrama- 



nian 



291 ]. which we assume here. These equations for Ml and H, describing the evolution 
of magnetic correlations at small and large scales. The effective diffusion coefficient k,n 
includes microscopic diffusion (rj), a scale-dependent turbulent diffusion [T L l (0) — T L l (?")] , 
and a ambipolar drift 2aM^ (0, t) , which is proportional to the energy density of the fluctu- 
ating fields. Similarly, aN is a scale-dependent a effect, proportional to [2C(0) — 2C (r)]. 
The nonlinear decrement of the a effect due to ambipolar drift is 4aH (0,t) , proportional 



to the mean helicity of the magnetic fluctuations. The G (r) term in equation a 



rapid generation of small scale magnetic fluctuations due to velocity shear 29 



30 



lows for 



37|, 



This turbulent spectrum simulates Kolmogorov turbulence 39j. As in the galactic inter- 
stellar medium, the protogalactic plasma is expected to have Kolmogorov-turbulence, driven 
by the shock waves originating from the instabilities, associated with gravitational collapse. 

In the galactic context, we can neglect the coupling term a^H as a very good approxi- 
mation since it is very small and consider only the evolution of Ml 29]. 

For turbulent motions on a scale L and a velocity scale v, the magnetic Reynolds number 
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(MRN) is R m = vL/t]. There is a critical MRN, R c « 60, so that for R m > R c [29(, modes 
of the small scale dynamo can be excited. The fluctuating field, correlated on a scale L, 
grows exponentially with a growth rate T L ~ u/L {^j]. 



IV. THE PARAMETERS OF THE TURBULENT PLASMA 



We use the fiducial parameters, suggested in the literature for the plasma that was present 
in the protogalaxy 40|, |4l|: total mass M ~ 10 12 M Q , temperature T ~ 10 6 K, and size 
L c ~ 200 kpc. The ion kinematic viscosity is ~ 5 x 10 26 cm 2 /s, the Spitzer resistivity r) s = 
6.53 x 10 12 T~ 3 / 2 In A cm 2 s~ 1 ~8x 10 4 cm 2 s~ 1 , and the typical eddy velocity V c ~ 10 7 cm/s. 



V. RESULTS 



In Fig. 1, we evaluate Ml for various values of r and in Fig. 2 for various values of 
V c , solving numerically equation (JSJ). In Fig 3 we evaluate the mean value of the magnetic 
field as a function of r and t. Our previous work jsj showed that the natural fluctuations of 
the primordial plasma predicted by the Fluctuation-Dissipation Theorem produces a cosmic 
web of randomly oriented dipole magnetic fields. The average field over a region ~ 2 kpc 
is predicted to be 0.003 nG. We assume this seed field and examine its amplification by 
the small scale dynamo in s protogalaxy. This seed field corresponds to an M^(~ B 2 ) ~ 
10~ 23 G 2 . Of particular interest is thus the growth of Ml with an initial value Mlo ~ 
10~ 23 G 2 in Figs. 1 and 2, for initial magnetic fields B ~3x 10~ 12 (2/cpc/r) 3/2 G of size r 
in Fig. 3. 



VI. CONCLUSIONS AND DISCUSSION 



It was shown previously that the magnetic fields, created immediately after the quark- 
ladron transition, produce relatively intense magnetic dipole fields on small scales at z ~ 10 
^j. We show here that the predicted seed fields of size ~ 2 kpc and intensity 0.003 nG at 
z ~ 10 can be amplified by a small scale dynamo in protogalaxies to intensities close to 
observed values. In the small scale dynamo studied, we use the turbulent spectrum given 
by Subramanian j^j. The characteristic velocity V c and length L c , used in the expression 
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t (years) 



FIG. 1: Values of Mx(G 2 ) as a function of t (years) and r. Solid black line has the reference values: 
M L (r,0) = 10- n (0.1pc/ r ) 3 G 2 , L c = 200 kpc, r = 3 kpc, and V c = 10 7 cm/s in Eqs. (36)-(38). 
Dashed red line is for r = 4 kpc, and the dotted blue line for r = 5kpc. 



for the vorticity V c /L c , are V c 
to that found by Ryu et al. 



7 cm/s and L c ~ 200 kpc. This vorticity is comparable 
, studying the formation of large scale structures. The 



length L c ~ 200 kpc used is a characteristic size o 



a protogalactic cloud. The turbulent 



spectrum used simulates Kolmogorov turbulence 39[. From our Figs. 1 and 2, we find that 
Ml(~ B 2 ) increases from ~ 10~ 23 G 2 (corresponding to a magnetic field B ~ 3 x 10~ 12 G 
over a region L ~ 2 kpc) to Mx ~ 10 18 G 2 (corresponding to a field ~ 10~ 9 G over a 
region L ~ 2 kpc) in 10 9 years. This corresponds to a ~ 6 e-fold amplification of B in a 
relatively short time. Collapsing to form galaxies at redshift z ~ 10, the density increases 
by a factor of ~ 200 and the magnetic fields are amplified by a factor of ~ 34. This predicts 
0.03 fiG fields over 0.34 kpc regions in galaxies. If the dipole magnetic fields predicted by the 
Fluctuation-Dissipation Theorem are not completely random, microgauss fields over regions 
> 0.34 kpc are easily obtained. The model studied is thus a strong candidate to explain the 
[iG fields observed in high redshift galaxies. 
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FIG. 2: Values of Mi(G 2 ) as a function of t (years), varying V c . Solid black line has the reference 
values in Fig. [TJ Dashed red line is for V c = 8 x 10 6 cm/ s. Dotted blue line is for V c = 6 x 10 6 cm/s. 
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FIG. 3: Values of B(G) as a function of t (years) and r(kpc) for reference values of Fig. 1. 
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